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Abstract 
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1 Introduction 


The interplay between a group action and the spectral analysis of a 
differential operator invariant under the action is a popular theme in 
analysis. If the group acts on a metric graph, the operators —D'^ + q 
with invariant q are obvious candidates for a spectral theoretic anal¬ 
ysis. This work treats operators —D^ + q on & metric Cayley graph 
Tm of the nonabelian free group Fm on M generators. These Cayley 
graphs are regular trees, with each vertex having degree 2M. In the 
present work the M edge types of Tm associated to the generators 
of Fm may have different lengths, with even functions q varying over 
the M edge types. Remarkably, novel techniques show that there is 
a system of M multipliers /im(A), resembling those of Hill’s equation 
[T6] , which can be used to construct the resolvents of the operators. 
Echoing the Hill’s equation analysis, the location of the spectrum is 
encoded in the behavior of the multipliers on the real axis. 

There is a large literature treating various aspects of analysis on 
symmetric inhnite graphs. Homogeneous trees were considered as dis¬ 
crete graphs in |3]. The quantum graph spectral theory of —D^ + q 
on homogenous trees was studied in [2], assuming that each edge had 
length 1, and that q{x) was the same even function on each edge. 
These assumptions meant that the graph admitted radial functions, a 
structure which facilitated a Hill’s equation type analysis of the spec¬ 
tral theory. The spectral theory of radial tree graphs was considered 
in [B] and [12]. A sampling of work exploiting this structure includes 
[31 El El El ED . Certain physical models can also lead to graphs of lat¬ 
tices in Euclidean space where the group (e.g. is abelian [igiiH]. 

The quantum Cayley graph analysis begins in the second section 
with a review of quantum graphs and the dehnition of the self-adjoint 
Hilbert space operator A + q which acts by sending / in its domain to 
— D’^f+qf. The third section reviews basic material on Cayley graphs, 
and in particular the Cayley graphs Tm of the free groups Fm- For each 
edge e of the Cayley graph and each A G C \ [0, oo), a combination of 
operator theoretic and differential equations arguments identifies a one 
dimensional spaces of ’exponential type’ functions which are initially 
dehned on half of Tm- The translational action of generators of Fm 
on subtrees of Tm induces linear maps on the one-dimensional spaces 
of exponential functions, thus producing multipliers fJ^mW for m = 
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1,..., M. A square integrability condition shows that |/im(A)| < 1 for 
all A G C \ [0, cxd). 

The fourth section starts by linking the multipliers and rather ex¬ 
plicit formulas for the resolvent of A -|- g. Recall that the multipliers 
for the classical Hill’s equation satisfy quadratic polynomial equations 
with coefficients which are entire functions of the spectral parame¬ 
ter A. In this work the multipliers /im(A) satisfy a coupled system 
of quadratic equations with coefficients that are entire functions of A. 
An elimination procedure shows that the equations can be decoupled, 
leading to higher order polynomial equations with entire coefficients 
for individual multipliers /im(A). The multipliers /Xm(A) have exten¬ 
sions from above and below to real a. The extension is generally holo- 
morphic, but as in the classical Hill’s equation the difference Smic) of 
the limits from above and below can be nonzero. Except for a discrete 
set, the spectrum of A -|-g is characterized by the condition Sm(o') ^ 0 
for some m. 

In the hnal section the system of multiplier equations is explicitly 
decoupled for the case M = 2. Computer based calculations are used 
to generate several spectral plots. 

2 Quantum graphs 

Suppose T is a locally hnite graph with a hnite or countably inhnite 
vertex set V and an directed edge set £. In the usual manner of 
metric graph construction [2], a collection of intervals {[0,/e],e G £} 
is indexed by the graph edges. Consistent with the directions of the 
graph edges {v,w), the initial endpoint v is associated with 0, and 
w is associated with /g- Assume that each unordered pair of distinct 
vertices is joined by at most one edge. As a result, the map from 
the directed graph to the undirected graph which simply replaces a 
directed edge (u, v) with an undirected edge [u, n] is one-to-one on the 
edges. A topological graph results from the identihcation of interval 
endpoints associated to a common vertex. 

The Euclidean metric on the intervals is extended to a metric on 
this topological graph by dehning the length of a path joining two 
points to be the sum of its (partial) edge lengths. The (geodesic) 
distance between two points is the inhmum of the lengths of the paths 
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joining them. The resulting metric graph will also be denoted T. 

To extend the topological graph T to a quantum graph, function 
spaces and differential operators are included. A function / : T —)■ C 
has restrictions to components fe ■ [0,/e] —t C. Let T^(r) denote the 
complex Hilbert space ©eL^[0,/e] with the inner product 

{f,9) = [ f9 = Y] [ fe{x)ge{x) dx. 

Jr , Jo 

Given a bounded real-valued function g on T, measurable on each 
edge, differential operators —D'^ + q are dehned to act component¬ 
wise on functions / G L^(r) in their domains. The functions q are also 
assumed to be even on each edge, qe{le —x) = qe{x). This assumption 
plays an important role as the analysis becomes more detailed. 

Self-adjoint operators acting by —D'^ + q can be dehned using stan¬ 
dard vertex conditions. The construction of the operator begins with a 
domain Vcom of compactly supported continuous functions / G T^(r) 
such that /' is absolutely continuous on each edge e, and f” G T^[0, lg\. 
In addition, functions in Vcom are required to be continuous at the 
graph vertices, and to satisfy the derivative condition 

5^a,/e(t;)=0, (2.1) 

where e ~ n means the edge e is incident on the vertex n, and = 
d/dx in outward pointing local coordinates. 

Since the addition of a constant will make the potential nonegative, 
but have only a trivial effect on the spectral theory, the assumption 

g > 0 (2.2) 

is made for convenience. With the domain T>comi the operators —D'^+q 
are symmetric and bounded below, with quadratic form 

{(-D^ + q)fJ)= j^\ff + q\f\\ (2.3) 

These operators always have a self-adjoint Friedrich’s extension, de¬ 

noted A -I- g, whose spectrum is a subset of [0,cxd). When the edge 
lengths of T have a positive lower bound the Friedrich’s extension is 
the unique self adjoint extension la p. 30]. 
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Say that an edge e = (n_,n+) G of a connected graph F is a 
bridge if the removal of (the interior of ) e separates the graph into 
two disjoint connected snbgraphs. If e is a bridge, let F^ denote the 
closnre of the connected component of F\nzp which contains the vertex 
v±. Less formally, F^ inclndes e, the vertices and the side of F. 

For A G C \ [0, cxd), the resolvents R{\) = (A + g — \I)~^ of the 
self adjoint operators A + g provide special solntions of —y” + qy = \y 
on F^. Let denote the space of fnnctions y± : Tf ^ C which (i) 
satisfy 

- y" + qy = ^y (2.4) 


on each edge e G F^, (ii) are continnons and sqnare integrable on F^, 
and (iii) which satisfy the derivative conditions (2.1) at each vertex of 
Fj^ except possibly n=p. 


Lemma 2.1. Suppose e is a bridge of the connected graph F, and 
A G C \ [0, cxd). Then Xf is one dimensional. 

Proof. Snppose two linearly independent fnnctions gi, g 2 on F+ satisfy 
(i) - (iii). Then a nontrivial linear combination h = aigi + a 2 g 2 
wonld satisfy h(n_) = 0. Consider altering the domain of A + g by 
replacing the vertex conditions at n_ by the condition /(n-) = 0 for 
each edge incident on n_. The resnlting operator is still self-adjoint 
and nonnegative on L^(F), and restricts to a self-adjoint operator on 
L^(Fj“). The fnnction h is then a sqnare integrable eigenfnnction with 
eigenvalne A, which is impossible. A similar argnment applies to F“. 
Thns X^ is at most one dimensional. 

Snppose ^ is a nontrivial solntion of the eqnation 


— z” ri qz = \z 


(2.5) 


on the interval [0, and / is a fnnction in the domain of A -|- g with 
snpport in e. As a fnnction on [0, lf\ the fnnction / then satishes 
/(O) = fife) = 0 and /'(O) = f\le) = 0. Integration by parts and the 
vanishing bonndary conditions for / lead to 

0= f f[—z" + qz — \z]dx= f [—f” + qf — \f]zdx. 

Jo Jo 

Since A is not in the spectrnm of A -|- g, the resolvent R{X) = 
{A + q — XI)~^ maps L^(F) onto the domain of A -|- g. Extend the 
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functions z by zero to the other edges of F to obtain an element of 
L^(r). Suppose / = R{X)z G L‘^{T) had support in e; then the above 
calculation would give 

pie 

/ dx = 0, 

Jo 

which is impossible. Also, since does vanish outside of e, the function 
/ = R{X)z satishes —/" + qf = Xf on each edge other than e. 

For i = 1,2 there are two independent solutions Zi of (2.5) on e, and 
after extension of Zi by zero, there are two independent functions fi = 
R{X)zi. As noted above, the functions ft must be nonzero somewhere 
on r \ e; without loss of generality suppose /i is not identically zero 
on r+ \ e. 

An argument by contradiction shows that at least one of the func¬ 
tions /i, /2 must be nonzero F” \ e. Suppose both fi and /2 are iden¬ 
tically zero on F“ \e. Dehne functions Qi G X+ which agree with fi on 
F+ \e, but which satishes (2.4) on e, with initial data gi{v+) = /+(u+) 
and with d^gi{v+) chosen so the derivative conditions ( |2.1 ) at are 
satished for gi. Since X+ is at most one-dimensional, a nontrivial lin¬ 
ear combination g = a\gi + a 25'2 = 0 on F+, and so / = ai/i -|- 02/2 
is zero on F \ e. But the existence of a nontrivial function / = R{X)z 
with support in e was ruled out above. 

The space X+ is then the span of gi, and the case of X“ is similar. 

□ 


The construction of Lemma 2.1 also provides the next result. 


Lemma 2.2. Suppose e is a bridge of the connected graph F, and 
A G C \ [0, 00 ). A basis g{X) ofXf may be chosen holomorphically in 
an open disc centered at X, and real valued if X E (—oo,0). If x ETf 
then g{x, X) is holomorphic. 


Proof. For A G C \ [0, 00 ) the resolvent R{X) is a holomorphic operator 
valued function, so the functions / = R{X)z are holomorphic with 
values in the domain of A-|-g. For a: G F=*"\e the evaluations /(x), f\x) 
are continuous functionals na p. 191-194] on the domain of A -|- g, 
so the values g{v+) and dugi{v+) from Lemma |A 1 are holomorphic, as 


are the L^(F^) functions g{X) and the values g{x, A) for x G e. All of 
these functions can be chosen to be real valued if A G (— cxd, 0). □ 
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3 Quantum Cayley graphs 

3.1 General remarks 

Suppose G is a finitely generated group with identity i. Let S C G be 
a finite generating set for G, meaning that every element of G can be 
expressed as a product of elements of S and their inverses. Following 
[T7] . the Cayley graph F^^g for the group G with generating set S is 
the directed graph whose vertex set V is the set of elements of G. The 
edge set of Fc^g is the set £ of directed edges (n, vs) with s G §, initial 
vertex n G G and terminal vertex vs. When confusion is unlikely we 
will simply write F for F^ g. Assume that if s G §, then ^ §. This 
condition avoids loops {v,vl), and insures that at most one directed 
edge connects any (unordered) pair of vertices. We will often consider 
F to have undirected edges [u,us], with the directions given above 
available when needed. 

G acts transitively by left multiplication on the vertices of F; that 
is, for every u, tn G V there is a G G such that w = gv. If e = 
{v,vs) G £, then ge = {gv,gvs) E £, so G also acts on £, although 
this action is not generally transitive. Say that two directed edges 
61,62 are equivalent if there is a G G such that 62 = 5 ' 6 i. The 
equivalence classes will be called edge orbits of the G action on £. 

Proposition 3.1. A set B of directed edges is an edge orbit if and 
only if there is a unigue s G § such that B = {{v,vs),v G G}. It 
follows that the number of edge orbits is the cardinality o/S. 

Proof. If u, tn G G, then wv~^{v, vs) = {w, ws), so for a fixed s G § all 
edges of the form (u, vs) are in the same orbit. If g(v, vsi) = {w, WS 2 ), 
then gv = w and gvsi = WS 2 , so wsi = WS 2 and si = S 2 . □ 

Proposition 3.2. IfG is generated by the finite set S, then the undi¬ 
rected graph F(g_s is path connected. 

Proof. If our requirements on generating sets are momentarily relaxed 
and S is extended to the set Si = {s, s G S}, then the Cayley graph 
Fij gi will have a directed path from every element of G to i. An edge 
of this graph has one of the forms (u, vs) or (u, As an undirected 

edge, [u,us“^] = [us“^,u] = [us“^,so for every directed edge 
of Fg,Si there is an undirected edge of Fg,s with the same vertices. 
Consequently, the undirected graph Fq^s is path connected. □ 


7 


Cayley graphs T can be linked with the spectral theory of differen¬ 
tial operators. To maintain a strong connection with the gronp G, the 
edges of T in the same G orbit will have the same length. The action 
of G on the combinatorial edges may then be extended to the edges 
of the metric graph by taking x G {0,l(^y^ys)) to gx = x E (0, 

This gronp action also provides a G action on the fnnctions / on T. 
The action simply moves the edge index, so that in terms of fnnction 
components gfge{gx) = gfge{x) = feix). Fnnctions are G-invariant if 
fge{x) = fe{x) for all directed edges e and all g E G. A qnantnm Cay¬ 
ley graph will be a qnantnm graph whose nnderlying combinatorial 
graph is the Cayley graph of a hnitely generated gronp, whose edge 
lengths are constant on edge orbits, and whose differential operator 
A -|- g commntes with the gronp action on fnnctions. Since there is 
little chance of confnsion, the same notation, e.g. T, will be nsed for 
the corresponding qnantnm, metric, and combinatorial graphs. 

3.2 Free groups and their graphs Tm 

This work will focns on Cayley graphs with G = Fm, the free gronp 
ini with rank M. Recall that the elements of are eqnivalence 
classes of hnite length words generated by M distinct symbols Si,... ,sm 
and their inverse symbols ; ■Sm • Two words are eqnivalent if 

they have a common rednction achieved by removing adjacent symbol 
pairs SmS^ or The gronp identity is the empty word class, the 

gronp prodnct of words Wi, W 2 is the class of the concatenation W 1 W 2 , 
and inverses are formed by nsing inverse symbols in reverse order, e.g. 

= Sh^S3S2^ 

Given a free gronp ¥m with generating set § = {si,..., sm}, Rt T/w 
denote the corresponding Cayley graph. These (nndirected) graphs 
(see Figure have a simple structure [I?! p. 56]. 

Proposition 3.3. The undirected graph Tm is a tree whose vertices 
have degree 2M. 

Proof. Suppose Tm had a cycle with distinct vertices tci,..., wk and 
edges [wk^wi] and [wk,Wk+i\ for A; = 1,..., iF — 1. In the undi¬ 
rected graph Tm edges extend from Wk by some Sm or so that 
Wk+i = WkSra or Wk +1 = WkS^, SO each vertex has degree 2M. The 
sequence of visited vertices Wi,... ,wk,Wi is described by a word of 
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Figure 1: A rescaled graph T 2 


right multiplications by the generators and their inverse symbols equal 
to i in ¥m- Since this word can be reduced to the empty word, it 
must have adjacent symbols SmS^ or This means the ver¬ 




tices tCi,..., wk are not distinct, so no such cycle exists. Since 7 m is 
connected by Proposition]^ and has no cycles, 7 m is a tree. □ 


3.3 Abelian subgroups and multipliers for Tm 

Each edge e = {v,vs) G Tm is a bridge. With u = u_ and vs = u_|_, 
the subgraphs P^ described above will be subtrees of Tm, denoted by 
T^. The vector spaces are as in Lemma 

Lemma 3.4. Suppose e = {v,vs) is an edge oJTm and X G C\[0, cxd). 
If is a nontrivial element ofTf, then is nowhere vanishing on 

T^. 

Proof. Suppose |/^(a;o) = 0 for some Xq G T^. First notice that y^ 
must then vanish identically on the subtree To consisting of points 
xi of with the property that paths from xi to v must include 


2.1 
























Xq. Otherwise, a nonnegative self-adjoint operator A + q could be 
obtained on by using the boundary condition /(xq) = 0. This 

operator would have a nontrivial square integrable eigenfunction, the 
restriction of to To, with the eigenvalue A G C \ [0,cx)), which is 
impossible. 

Since solutions of —y" + qy = \y are identically zero on an edge e 
if y^xo) = y'^xo) = 0 for some xq G e, we may assume xq is a vertex. 
Since the function y^ vanishes identically on To, the continuity and 
derivative conditions at xq force y to vanish on all the edges with xq 
as a vertex. The function y^ must now be identically zero on T^^, 
contradicting the assumption that the function was nontrivial. □ 

The structure of the elements of is strongly constrained by the 
symmetries of combined with the fact that X^ is one dimensional. 
A simple observation is the following. 

Lemma 3.5. Suppose ei = {vi^viSm), S' G G, and 62 = (v2,V2Sm) = 
(gvi,gviSm)- Forj = 1,2 let yj G X^^. withyj{vj) = 1. Then gyi = y2. 

Proof. The action by g is an isomorphism of and Since X^^ 
is one dimensional, gyi is a scalar multiple of 1 / 2 • These two functions 
agree at V 2 , so are equal. □ 

For each vertex v and integers k, left multiplication by the abelian 
subgroup of elements acts on Tm- These maps carry the edge 

e(0) = (n, vSm) to the edges e{k) = vs^^). The key role of these 

group actions is related to the following geometric observation. 

Lemma 3.6. The trees T^|^^ are nested, with C T^p.)- In addi¬ 

tion, 

00 

I'm =(J I'e(-k)- 

k=0 

Proof. Other than vs^, the vertices of the trees are those el¬ 
ements of Fm which have a representation vs^Smg, where s^g is a 
reduced word in Fm- If w = vs^^Smg G is a vertex with Smg 

reduced, then s'^g is reduced and w = vs^SmSmg G le^k)- Thus the 
trees are nested. 

More generally, for any integer k, a word w & may be repre¬ 
sented as w = vs^Smg with Smg = si^v~^w. First take a reduced 
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representative u for v~^w. Suppose u begins on the left with a string 
followed by an element of S different from Taking k = j — 1 
gives the desired form, and every vertex w is in some 


□ 


Suppose y G X^(o) satishes y{v) = 1, and 2 ; G satishes 

z{vSm) = 1 - Since C 'T'^oy the restriction of y to is an 
element of Because y is nonvanishing, there is a nonzero multi¬ 
plier G C associated to each generator Sm such that y = 

on Tyy. In particular /im(A) = y{lm)/y{0). 

Lemma 3.7. The multipliers fimW o,re holomorphic for A G C \ 
[0, cx)), with Pm(A) = /im(A). 


Proof. By Lemma 2.2 the formula /i. 
(A) is holomorphic when A G 


pi. 


(A) = y{lm)/y{Q) 
C \ [0,oo). If A G (■ 


shows that 
- 00 , 0) and 


y is chosen real, then /i(A) is real. 


The two functions pim(^) and 
(A) are holomorphic and agree for A G (—cxd,0), so agree for all 

□ 


/^m 

A G C \ [ 0 , 00 


Because the function q is even on each edge, that is q{lm — x) = 
q{x), the same multipliers will arise when comparing elements of X“ 
if the edge directions are reversed by using the generators of 
instead of Sm- These multipliers provide a global extension of functions 

ill ^5o)- 

Lemma 3.8. Suppose C C If yj G with yj{v) = 
1, and yk G ’^^k) 2/fc('^) = 1? then yj = yk on %(j). Elements y"^ 

o/X^^Qj may be extended via the multipliers to functions defined on all 
ofTM- 


Proof. The function yk restricts to an element of Since nontrivial 
elements of never vanish, but yk{v) — yj{v) = 0, the difference is 
the zero element of Since these extensions of X^^q^ are consistent, 
elements y^ of X^^^^ extend via the multipliers to functions dehned on 
all of the trees Tyy. □ 

Lemma 3.9. For A G C \ [0, 00 ), the multipliers satisfy |/im(A)| < 1. 
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Proof. Recall that y is nowhere vanishing, so 


blVo. 


A nontrivial element y of is square integrable on so in partic¬ 
ular 


>+i 


5^/ \y\^= 


k=0 


k=0 


and |/r^(A)| < 1. 


□ 


4 Analysis of the multipliers 

On each edge [0, Im] the space of solutions to the eigenvalue equation 


(2.4) has a basis Cm{x, A), Smix, A) satisfying 0^(0, A) = 1 = 5^(0, A) 


and 5^(0, A) = 0 = C'^{0, A). These solutions satisfy the Wronskian 
identity 

Cm{x, X)S'Jx, A) - S'^ix, X)CUx, A) = 1. (4.1) 

If g = 0 and ui = a/A, these functions are simply Cm{x, X) = cos(a;x), 
Sm{x, A) = sin(a;x)/a;. 

If A) = 0 then A is an eigenvalue for a classical Sturm- 

Liouville problem, implying A G [0, cxd). For A G C \ [0, cxd) there 


is a unique solution of (2.4) with boundary values |/m(0. A) = a, 
ym{lm,X) = given by 

y^{x,X) = aCm{x,X) + (4.2) 

Because qm{x) = qm{}m — x) for each edge, there is an identity 
Cm(^m - X, A) = A)Om(x, A) - A)S'm(x, A) 


since both sides of the equation are solutions of (2.4) with the same 


initial data at x = Im- Setting x = 0 leads to the identity 

C'^(/^,A) = R:„(/^,A). (4.3) 


In addition to the coordinates originally given to the edges of 7m, it 
will be helpful to also consider local coordinates for 7^+ which identify 
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edges with the same intervals [0,/^], but with the local coordinate 
increasing with distance from a given vertex v. Since q is assumed 
even on each edge, the operators A + q are unchanged despite the 
coordinate change. 

4.1 Multipliers and the resolvent 

The next results show that edges in the same orbit have the same 
multipliers. 

Theorem 4.1. Assume e = {v,vsm), A G C \ [0, cx)), and y G X+ 
with y{v) = 1. Suppose the edge e in is in the same edge orbit as 
e, with the local coordinate for e increasing with the distance from v. 
Using the identifications of e and e with [0,/^], the restriction yi of y 
to e satisfies 

yiilm) _ y{lm) _ /.\ 

2/1 ( 0 ) 2 /( 0 ) 

Proof. If w is the vertex of the edge e closest to v (see Figure]^, then 
e has one of the forms {w,wsm) or (ws"^, w). In the hrst case, where 
wv~^e = e, the tree is a subtree of and translation by wv~^ 
carries to X+. As functions on [0,/^], 2/i is a nonzero multiple of 
y since X+ and X+ are one dimensional. 

In the second case, when e = {wsf^,w), the tree is generally 
not a subtree of 71'*“, but TP _i , is. A different argument will reduce 
the second case to the hrst. As undirected graphs there are isomor¬ 
phisms between the trees TA,,,, n and Tp _i .. One such is obtained 

f {w,WSm) (wsP,w) 

by interchanging the roles of Sm and There is a corresponding 
involution of X^ obtained by interchanging function values on the iso¬ 
morphic trees. Since X+ is one dimensional, this involution is given 
by a constant factor. The nonzero value of y at the vertex w is hxed 
by the involution, so the tree interchange must leave the functions 
hxed. □ 

Corollary 4.2. Assume e = {v, vSm), A G C\ [0, oo), and y G X+ with 
y{v) = 1. Suppose that for j = 1,2, the edges Cj = {wj^WjsP) in 
are in the same edge orbit, with the local coordinates for cj increasing 
with the distance from v. The restrictions yj of y to Cj satisfy 

yiilm) _ yiilra) _ /.\ 
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WSm 



V 


VSm 


Figure 2: Edges in a common orbit 


Proof. If Cl = u), then is a subtree of Tff. Since ei and lie 

in the same edge orbit, the previous theorem may now be applied. □ 


As a consequence of Theorem |4.1 and Corollary 4.2| the functions 
y G have the following description. 


Theorem 4.3. Assume A G C \ [0,cxd), e = {v,vSm) and y+ G X+ 
with y+{v) = 1. Suppose w is a vertex in , and the path from v to 
w is given by the reduced word . Then 


pmPk(l) ■ ■ ■ Pk(n)- 


(4.4) 


Using (4.2), the vertex values of y+ can be interpolated to the edges. 


Because the functions in X+ are continuous, the multipliers /ifc(A) 
are simply the value at of the solution y^ in X+ with initial value 1 
at X = 0 on edges of type k. That is. 


Pk{^) — Ck{lk, A) + |/fc(0)S'fc(4, A). 


(4.5) 


Theorem 4.3 may also be used to describe the functions ?/_ G Xg . 
The functions y+,y- can be used to construct the resolvent R{X) = 
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[A + g — A]“^ on C \ [0, oo). If the (nonvanishing) functions and ?/+ 
were linearly dependent on e, then there would is a nonzero constant 
c such that y-{x) = cy+{x) for a; G e, and the function 

(y-ix), xeT-,-) 

\cy+{x), X G T+J 


would be a square integrable eigenfunction for A+q. Consequently, the 
functions y_ and y+ must be linearly independent on e if A G C\ [0, oo). 

In particular for each A G C \ [0, oo) the Wronskian hFfc(A) = 
y+y'_ — y'+y- for edges of each type k is nonzero, and independent of 
X. By using (4.5) the Wronskian Idfc(A) can be expressed in terms of 
the multipliers. Consider evaluation of fhfc(A) at x = Ik- Compared 
to y+, which satishes (4.5), with g+(0. A) = 1, the function y_ would 
have the edge direction reversed. This function has y-{lk) = 1, and 
because of the reversed edge direction. 


y-{x, A) = Ck{lk -x,A) - y'_{lk)Sk{lk - x, A), 


so that 

hfc(-^) = 2/-(0, A) = Ck{lk, A) — yL{lk)Sk{h, A). 
Evaluation at x = Ik gives 


Wk{\) = {y+y'_ - y'+y-){lk) 


/\\^k{^ky^) k'ki^) \r^l n \\ I Ck{lki Q/ -.N-i 

- -FTTl--H-FTTl-- ^kil-k, AjJ 

and the identities ( 4.1[ ) and (4.3) give the simplihcation 

l-h^(A) 


Sk{lk,X)' 

For A G C \ [0, oo) and t G [0, Zg], dehne the kernel 
f> A)|/+(f, X)/Wk{X), 0<x <t<le, 

Ke[X,t,A) y_^t^X)yAx,\)IWt{>-), 0<t<X<k. 


(4.6) 


(4.7) 


The following observations show that the values of x can be extended 
to the whole of 7m- 
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If /e is supported in the interior of e the function 


he{x) = / Re{x,t,\)fe{t)dt 

Jo 

satishes [A + q — X\he = fe, and in neighborhoods of 0 and Ig the 


function hg satishes (2.4). The kernel Rg{x,t,X) and the function hg 
can then be extended to Tm using the values of and y+ on . The 
extended function hg is square integrable on Tm and satishes the vertex 
conditions, so hg agrees with the image of the resolvent acting on /e, 
that is hg = R{X)fg. Since the linear span of functions fg is dense in 
L‘^{Tm)i and the resolvent is a bounded operator for A G C \ [0,oo), 
the discussion above implies the next result. 

Theorem 4.4. For A G C \ [0, cx)), 

= r Re{xT\)fe{t)dh feL^TM), 


the sum converging in L‘^{Tm] 


4.2 Equations for the multipliers 

Theorem 4.5. For A G C \ [0, cxd) and m = 1,..., M , the multipliers 
Pm (A) satisfy the system of equations 

hliW - 1 _ 2 ^ Pfc(A) - Ckjlk, A) ^ . . 

>S'm(^m, A)pm(A) “ A) 

Proof. Begin with an edge e = (u,usm)- In addition to e, the ver¬ 
tex vSm has 2M — 1 other incident edges. One is the type m edge 
(uSm, while the others have one of the type k forms (uSm, vSmSk) 
or {vSmSf^,vsm): where k G \ {m}. Let y G X+ satisfy 

y{xsm) = 1 , and let yk, respectively ym denote the restriction of y to 
the subtrees with root vertex vSm and initial edges of type k, respec¬ 
tively m. As a consequence of Corollary |4.2| , for a hxed value of k the 
two functions yk agree as functions of the distance from vSm on the 
two edges of type k incident on vSm- 

Let Zm denote the restriction of y to the edge e. The continuity 
and derivative vertex conditions at vSm relate Zm to the restrictions 


16 





ykillm- Using local edge coordinates which identify vs^ with for 
the edge e, and which identify vsm with 0 for the other incident edges, 
the initial data for Zm at vSm is 

Zmilm, A) = 1, A) = A) + 2 ^ |/fc(0, A). 

k^m 

The function Zm{x, A) may be written as 

Zjn{,x, A) = - X, A) - [|/^(0, A) + 2 ^ |/fc(0, A)]S'm(/m - x, A). 

kj^m 

Evaluation at x = 0 gives 

Zm{0, A) = Cm{lm, A) — [|/m(0, A) + 2 y'ki^-, A)]5'm(^m, A). 

k^m 


4(0, A) = -C(4. A) + 14(0, A) + 2 5^ 4(0, A)|4(4, A). 

k^m 

As noted above, on [0, Im] the function Zm{Xj A) is a scalar multiple 
of ym{x, A), so that Zm{x, A)/4m(0, A) = ym{x, A). In particular. 


l/m(0,A) 


^m(O.A) 

'^r7i(0. A) 


-c;(4,A) + |g;.(0,A) + V]4(4,A) 

Cm{lm, A) - [y'mi^, A) + 2 Y.k^m Vki^^ A)]S'm(/m, A) 
This can be rewritten as 


Cm{im, A)|/m(0, A) + A) 


(4.9) 


14(0, A) + 2 4(0, A)] [4(4, A)4(0, A) + 4(4. A)]. 

k^m 


Using (4.5) to substitute for 1/^(0) in (4.9) gives 
. h'nA(A) U*m(^m) A) 


Cmiuxy- 


+ CL{lm,X) 


SmiUX) 

r/^m(A) C*m(^ A) , r, yk{x) — Ck{ik, X) 


(4.10) 


Smilm,X) 


+ 25 : 


kj^m 


Skih, A) 


17 









X [f^mW — Cm{lrm -^) + A)]. 

With the help of the identities (4.1) and ( |4.3 ), these equations can be 
rewritten as (4.8). 

□ 


The solutions /ii(A),... ,/iM(A) of (4.8) coming from the resolvent 
of A + g can be recognized by a square summability condition. 

Theorem 4.6. Suppose A G C \ [0,cxd), e = {v,vSm), and w denotes 


a vertex in . Assume that /ii(A),... ,/im(A) satisfy (4.8). Define a 


function y+ : Tfi C by taking y+{v) = 1, defining y+{w) by (4.4) 
and interpolating the vertex values of y+ to the edges using (0. 
The function y^ is sguare integrable on Tfi if and only if 




(4.11) 




If y- : Tg C is defined similarly, and (4.11) is satisfied, then 


the formula (4.7) gives the resolvent of A + q as in Theorem 4-4 


Proof. The solutions of (4.8) satisfy pirai^) 0. For j = 1,2, two 
edges {wj, WjSm) G Xt have vertex values satisfying y+{wjSm)/y+{wj) = 
so the interpolated edge values yj given by (|4.2) satisfy 


2/2 (0)|/1 (x) = ?/i (0)7/2 (x), 0<x<k 


As a result. 


“IS 

so by comparing edges of type m with one of the edges incident on 
vsm we see that y+ is square integrable on Tfi if and only if (4.11) 
holds. 

Running the argument of Theorem 4.5 in reverse shows that the 
continuity and derivative conditions (2.1) hold at the vertices of Tfi' 
except possibly at v. If (4.11) holds, then //+ G X+ and the claims 
about the resolvent formula follow. □ 


The equations (4.8) have implications for the decay of /im(A). 
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Theorem 4.7. For 0 < a < tt, let Q denote the set of X with \ arg(A) — 
7 r| < a and |A| > 1. For A G ^2, 


lim |/im(A)|e 

I A|^oo 




2M 


(4.12) 


Proof. Rewrite (4.8) as 


M 


VAFm(A,/}7i) (A,(A) 


^ /^fc(A) - C*fc(A, Zfc) 


k=l 


VxSk{x,ik) 


Take Q'(\/A) > 0 and recall that |/ifc(A)| < 1. The functions 
Ck{X,lk), Sk{X,lk) satisfy the estimates [201 P- 13] 


|C'fc(A,/fc) - cos(yA/fc)| < 


C 


exp 


(|o(^)IW, 


\St{\h)- 


sm{'/Xlk), , C 


^ I < P^exp(|9(yA)4), 


while 


cos(v^a:) = -e-'^"(l + sin(yAa:) = ^e-'^"(l - 


For A e n, taking |A| —)■ oo in (4.13) gives 


lim 


lAl^oo y/XSm{lm)tJ^m{X) 


= -2iM, 


which implies (4.12). 


□ 


For each A G C the equations (4.8) are a system of polynomial 
equations Fm(^i, • • •, ^m, A) = 0 for m = 1,..., M which are satished 
by the multipliers pi,..., hm- The independence of the equations and 
local structure of the solutions may be determined by computing the 
gradients VF^ with respect to ..., with A treated as a param¬ 
eter. Recall from Lemma 3.9 that the multipliers satisfy \nm{X)\ < 1 
if A G C \ [0, cxd). 
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Theorem 4.8. Suppose X e C, Sm{X,lm) 7^ 0, and + 1 7^ 

0 for m = 1,... ,M. Then the complex gradients VPm are linearly 
independent if 

m=l 

Proof. The relevant partial derivatives are 

dPm _ I . I . 

dU SUlm)kliX) ^ 

and for j ^ m 

dPra _ ^ 1 

That is, there is an m-independent vector function W such that 


yPm = v^ + w, w 


^ 1/5'i(/i) ^ 

\1/ Sm{Im) / 


with Vm having m-th component equal to 


dPm/dfm + 2 / Smilrn) 


1 . 1 




+ 1 ] 


and all other components zero. 

If the vectors 14 + hh are linearly dependent, then there are con¬ 
stants Ofc not all zero such that + hh) = 0. If Smilm) 7^ 0 

the component equations can be written as 




1 


+ 1 ] 


M 


k=l 


This linear system is 


diag[ 





= 2 

/I .. 

• 

/ttl\ 

\OiM j 


•• 

• V 

\OCM ) 
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If none of the have the value —1, this system says [ai,..., q;m] 
is an eigenvector with eigenvalue 1/2 for the matrix 






€L\ 

2 , 1 




Cf+1 




Vectors with ^ a/c = 0 are in the null space of this matrix, and the 
remaining eigenvalue is the trace, with eigenvector 


Thus the condition for dependent gradients is 

Sn 


■ gf 
Tf+i’ 


M 


E 


*2 

sm 


^ +1 

=1 ^ 


= 1 / 2 . 


□ 


By applying the inverse and implicit function theorems for holo- 
morphic functions m p. 18-19] we obtain the following corollary. 

Corollary 4.9. Suppose A G C, Smi^Jm) 7^ 0, and 

Sn 


,^m(A) -h 1 7^ 0, 


M 


= 1 ^ 


m 


^1/2, m = l. 


,M. 


Then the solutions of the system (4.8) are locally given in C™ x C by 


a holomorphic valued function of X. 

Theorem 4.10. There is a discrete set Zq C M and a positive in¬ 


teger N such that for all X G C \ Zq the eguations (4.8) satisfied by 
the multipliers pirn{X) have at most N solutions £i(A),... ,£m(A). For 
A G C \ Zq, the functions frm{X) are solutions of polynomial eguations 
Pmi^m) = 0 in the one variable £m of positive degree, with coefficients 
which are entire functions of X. 

Proof. The polynomial equation in the single variable £i(A) will be 
considered; the equations satished by the other functions /im(A) may 
be treated in the same manner. 


Notice that the system (4.8) has the form 


Fi{^i, X) = 5 f(£i,... ,£m, A),..., Tm(£i, A) = g{^i,... ,£m, A), 
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with 


(4.15) 




e(A) 




M 


■ ■ ■ ,'Cm, A) = 2 


ik{F) — Ck{ik) 

■ 


k=l 

Subtraction of successive equations eliminates the right hand sides 
from M — 1 equations, giving a system of M equations, the equations 
indexed by the value of j = 1 ,..., M, 

J = 1 A) = ... ,^M, A) 

j = 2 Fi(6,A)-F2(e2,A) = 0, 


j — M A) — Fm{^m, A) — 0 

For m > 1, the m-th equation can be written as 

e-i(A)-i 


= (: 


' “Sm—I'Cm—1 (A) 

or by using the quadratic formula. 


+ (4.16) 


2e„4(A) - 


e-i(A) - 1 
“Sm—1'em—1 ( A ) 


5™ = 


'S'm—l^m—1 (A) 


5^+4 

^ m I ^ 


1/2 


(4.17) 


The variables • • • ,^2 can be successively eliminated from the 
hrst equation. Starting with k = M and continuing up the list of 
indices, the first equation can be written as a polynomial equation 
for with coefficients which are polynomials in • • • 5 ^^-! and the 


entire functions 5^(4, A). Repeated use of the substitution (4.16), 


followed by clearing of the denominators, reduces the hrst equation to 
degree one in These substitutions result in an equivalent system 
of equations as long as A ^ Zq, where 




(4,18) 


Now solve for 2^k, subtract [iCm-i(A) — l]*S'm/[*S'm-i'Cm-i(A)], use the 


substitution (4.17), and square both sides. Since squaring is a two- 


to-one map, it will not change the dimension of the set of solutions. 
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After applying these substitutions, the variable has been eliminated, 
and after clearing the denominators, the modihed j = 1 equation is 
a polynomial in with entire coefficients. The sub¬ 

stitution process provides a common bound N for the degrees of the 
polynomials Pm- 

Suppose the hnal version of the hrst equation does not have pos¬ 
itive degree for Dehne Qm = Fm{^rn, A) - Fm+i{^m+i, A) for m = 
1,..., M —1. The system of equations for ^i,..., is then the system 

Qi = 0,..., Qm-1 = 0, 

where 

dQm _ 1 r 1 .1 dQm _ “1 r 1 

and all other partial derivatives are zero. Suppose 0, 

and ^m(A)^ -|- 1 7 ^ 0 for m = 1,... ,M. Then the M — 1 gradients 
VQm are linearly independent, so outside of a discrete set of A the 
functions ^ 1 ,... would be holomorphic functions of that 

is, the solution set would have dimension 2, contradicting Corollary |4.9| 
which showed the dimension is 1. 

□ 


4.3 Extension of multipliers to [0, 00 ) and the spec¬ 
trum 


The mapping z —)■ A given by 


A = 


1 + z 


z = 


x/A 


1 + ^ 


is a conformal map from the unit disc {|z| < 1} onto A G C \ [0, 00 ). 
By using this conformal map and Lemma 3.9 the functions prn{X{z)) 


may be considered as bounded holomorphic functions on the unit disc. 
Classical results in function theory m p. 38] insure that Pm{X{z)) 
has nontangential limits almost everywhere as a function of z, and so 
the limits 

p^{a) = \im Praia ± ie) (4.19) 

£->0+ 
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exist almost everywhere on [0,cx)). By (4.8) and (4.18), the values 
are bounded away from zero uniformly on compact subsets of 

C\Zo. 

Since the functions /im(A) satisfy the polynomials equations Pmi^m) - 
0, more information about p^{cr) is available. The equations = 

0 have entire coefficients and positive degree for A G C \ [0, cxd). Let 
Zm C C denote the discrete set where the leading coefficient vanishes. 
A contour integral computation which is a variant of the Argument 
Principle, [H p. 152] or problem 2 of [121 P- 1^4], shows that for 
A G C \ Zm the roots of PmW, in particular PmW, are holomorphic 
as long as the root is simple. For A G C \ Z^ the roots extend con¬ 
tinuously to A G [0, oo) even if the roots are not simple. The limiting 
values /i^(cr) need not agree; let 

6m{(T) = aG[ 0 ,oo). 


Proposition 4.11. ///i+(cr) = /i^(cT) fora G {a,(3)\Zm, then 
extends holomorphically across {a, (3). 


Proof. On any subinterval (ai,/9i) C {a, (3) where PmW extends con¬ 
tinuously to the common value p^{a), the extension is holomorphic 
by Morera’s Theorem [Tnl p. 121]. The points in the discrete set 
Zm n (a, (3) appear to be possible obstacles to the existence of a holo¬ 
morphic extension, but since the extended function /im(A) is bounded 
the extension can be continued holomorphically across Z^C) {a, (3) too 
by Riemann’s Theorem on removable singularities. □ 


Theorem 4.12. Assume a ^ Zq. Form = 1,..., M suppose p^{a) ^ 
±1 and PmW extends holomorphically (resp. continuously) to a eM. 
from above (resp. below). Then the kernel function Re{x,t,3^) of ( 4.7| ) 
extends holomorphically (resp. continuously) from above (resp. below) 
to 


R: 


'di 


(x,t,a) = lim RJx,t,a + ie), 

e-S>0+ 


a G 


Proof. The Wronskian formula (4.6) shows that l/lTm(A) extends 
holomorphically (resp. continuously) if /im(A) does and 


Theorem 4.3 shows that the vertex values y±{w) extend in the same 


fashion as the multipliers pm- Finally, the interpolation formula (4.2) 


provides a holomorphic extension of y± from the vertex values as long 
as (T ^ Zo, that is Sm{lm, a) 7 ^ 0 for m = 1,..., M. □ 
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Recall [221 p. 237,264] that if P denotes the family of spectral 
projections for a self adjoint operator, in this case A + g, then for any 

/ e L^T) 


^[Pla,b] + Pia,b)]f = lim ^ ^ [R{a + ie) - R{a - ie)]f da. (4.20) 

Theorem 4.13. Suppose {a, (3) fl Zq = 0. Form = also 

assume that {a, (3) fl = 0 and that fi^{a) ^ ±1 for all a G (a,/9). 
If [a, b] C {a, (3), e is an edge of type m, and f G lP‘{e), then 


P\aMf = ^ y - Re{^)]f da. (4.21) 

If ai G (a,/d), then ai is not an eigenvalue of A + q. 

Proof. As noted above, the assumption that (a,/3) fl = 0 means 
the multipliers pm extend continuously to [a, h] from above and below. 
Since Pm(o') ^ ±1 the function extends continuously to [a, h]. 


Based on Theorem 4.3 and the interpolation formula (4.2), the kernel 
/?e(x, t, A) described in (4.7) extends continuously to [a, b] from above 


and below. The convergence of Re{x, t, a Pie) to Re{x, t, a) is uniform 
for t, X coming from a hnite set of edges. 

If the support of g G L'^{Tm) is contained in a hnite set of edges. 


then (4.20) and the uniform convergence of Reix, f, aPie) to Re{x, t, a) 
gives 


{\[P[aM P Pia,b)]f,g) = ^ ^ [Rei(^) “ (t^)](/,^) da. 

The set of g with with support in a hnite set of edges is dense in 
L‘^{Tm)i so the restriction on the support of g may be dropped. Sup¬ 
pose g is an eigenfunction with eigenvalue ai G (a, b) and with ||g|| = 1 , 
while / is the restriction of g to the edge e. Then the continuity of 
Re(x, t, a) means there is a Cg such that 

iRt(^) - Re((^)](f,ff) da\ < Ce\b-a\. 



This implies {P„^g,g) = 0, so the eigenfunction g doesn’t exist. Fi¬ 
nally, the absence of point spectrum in (a,/9) means that P[a,h\ = 
P{a,b), giving the formula (4.21). 

□ 
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Theorem 4.14. Assume a G [0, oo) \ Zq and for m = 1,..., M sup¬ 
pose Then a is in the resolvent set of A -\- q if and only 

if dm(o') = 0 in open neighborhood of a for m = 1,..., M. 


Proof. If a is in the resolvent set then the kernels described in (4.7) 


will have a common holomorphic extension to a from above and below. 
Evalnation gives 




l]Sm{lmA) 


SO that (/i+)^(cT) = (/i“)^(cT). A second evalnation, 
i?„(0,0,A) = 

1 -/^m(A) 

shows /i+( ct) - /U” (cr) = 5^{a) = 0. 

Now snppose Sm(o') = 0 in open neighborhood of a. Theorem 4.12 
notes that Rm{x,t,X) extends holomophically to a neighborhood of 
a. li g E L^iT) with snpport in a finite set of edges, the fnnction 
{Rm{X)f,g) also extends holomorphically as a single valned fnnction 
in an interval {a, ft) containing a. If [a,b] C {a, ft), then for any 
/ e L2(e) 

+ T{a,b)]f,g) = 0 . 

The set of g with with snpport in a finite set of edges is dense in T^(r), 
so l[P[a,b]+Ria,b)]f = 0 for all / e L'^ie). By linearity [P[a,6]+ -P(a,fe)]h = 
0 for any h G T^(r) with snpport in a finite set of edges, and since the 
projections are bonnded we conclnde that P[a,b\ + P{a,b) = 0 and (a, ft) 
is in the resolvent set da p. 357]. 

□ 


Corollary 4.15. Assume a G [0, cxd) \ Zq and for m = 1,...,M 
suppose p-m(o') 7^ ±1. Then a is in the resolvent set of A+q if and only 
if /im(cr) is real valued in open neighborhood of a for m = 1,..., M. 


Proof. If Pm(o') is real valned, then the symmetry pim{X) = /im(A) 
established in Lemma 3.7 means fif^icr) has the same real valne. The 
same symmetry also implies that Smicr) 7^ 0 if valned. 

□ 
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Theorem 4.16. For q > 0 and M > 2 the spectrum of A + q has a 
strictly positive lower bound. 

Proof. Since 


((A + g)/,/) = (A/,/)+ / g|/|2, 

JTm 

it suffices to verify the result when g = 0. 

Consider the case when the edge lengths Im are all equal to 1. Then 
the system (|4.8) reduces to 


(2M - l)/r^(A) - 2MC'(1, A)/i(A) + 1 = 0. 
The quadratic formula gives 

2MC(1, A) ± ^4M2C'2(1,A) -4(2M- 1) 


/i(A) = 


2(2M- 1) 


Since (7(1,0) = cos(O) = 1, the discriminant has the positive value 
4M^ — 4(2M — 1) = 4(M — 1)^ when A = 0, and p(ct) is real as long 
as cos^(A/a) > (2M — 1)/M^. 

Returning to the general case of a graph T with unconstrained edge 


lengths, recall (2.3) that the quadratic form for A is 


(A/,/)= / U'?- 

JTm 

Let a; be a coordinate for intervals [0,/m] and t for the interval [0,1]. 
For m = 1,... ,M let a; = 4>m{t) be a smooth change of variables. 
Assume (f' > Ci > 0 and 4>'{t) = 1 for t in neighborhoods of 0 and 1. 
If / is in the domain of A for the graph T, then / o 0 will be in the 
domain of A for a graph T i whose edge lengths are all 1. 

The chain rule and the change of variables formula for integrals 
give 


f*^m f* 1 

|2 


|/(a:)| dx= \f{(j)^{t))\ (j)'^{t) dt, 


and 


.df{x) 


dx 




0m (^) 
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As a consequence there is a constant C > 0 such that 


The calculation for graphs with edge lengths 1 shows that the expres¬ 
sion on the right has a strictly positive lower bound. □ 


Corollary 4.17. Suppose M > 2 and the lengths Im are rational. 
Then the resolvent set of A includes an unbounded subset of [0, oo). 

Proof. Assume A G C \ (— cxd, 0) so that y/X may be taken to be con¬ 
tinuous and positive for A > 0. In case q = 0, 

A) COs(/yj.2^^^): A) siu^/^V^), 


and these functions are periodic in a/A with period 27r/Im- If Im = 
Tm/rim with Tm, Vm positive integers, then the functions have a common 
period p = 27rnm=ihm- 

The functions C{lm,X) and S{lm,X) appear as coefficients in the 
equations (4.8). After multiplication by 1/\/A, the equations (4.8) 
exhibit the same periodicity, so have identical solutions for A and Ai 
whenever y/Xf = y/X -|- 2np for any positive integer n. 

By Theorem |4.6 the solutions of (4.8) which are multipliers are 
determined by the summability condition if A G C\[0, oo), so pm([V^+ 
2np]^) = Pm(A) for nonreal A. This identity extends by continuity to 
A G [0, oo). By Theorem 4.16 there is a cxo > 0 such that [0, cxo) is 
in the resolvent set of A. Corollary 4.15 shows that except possibly 
at a discrete set of points, the points a G [0, oo) which are in the 
resolvent set are characterized by real values of the multipliers 
so except for a discrete set of possible exceptions, (4n^p^, [^/(TQ + 2npY) 
is a subset of the resolvent set for A. □ 


5 Sample computations 


In this section some sample spectral computations are carried out for 
the case M = 2. The hrst step is to reduce the system of equations 
(4.8) to equations for individual multipliers. Two equations of degree 
four with entire coefficients are obtained. For g = 0 these equations are 
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solved numerically (using Matlab) for positive values a of the spectral 
parameter. After eliminating spurious solutions, the multiplier data 
is displayed in several hgures. 


5.1 Elimination step 


When M 


2 the system of equations (4.8) may be written as 


- 1 /ri(A) - C^ih, A) f,2{X) - C^ih, A) 
5i(/i,A)/xi(A) 5i(/i,A) ^ S.ikA) ’ ^ ^ 


^^l{X) - 1 /xi(A) - C^ih, A) fi2{X) - A) 
S2{h, X)fl2{X) Si{li,X) S2{l2,X) 

Subtracting the second equation from the hrst gives 

- 1 _ h2(^) - 1 ^ Q 

^l(/l,A)/ii(A) S2{l2,X)^I2{X) ’ 


Solving this quadratic equation for 2/i2(A) gives 


2/i2(A) — — 7 ^ 82 ( 12 , X) 


^i(h,A)/xi(A) 


(5.2) 


= ± 


(h^(A) - 1)^ 
ls!{h,x)^,l{x) 


A) + 4 


1/2 


(5.1) is already hrst order in /i 2 (A), and may be rewritten as 


= -2S,(k, + 2 C 2 (h, A). 

Oipi, A) 


Replacing the left hand side using (5.2) and squaring gives 

(h?(A) - 1)^ 


Sf{h,X)ftl{X) 


^2l/2,A)+4 = 452^(/2,A) 


/ii(A)-C'i(/i,A) 


5i(/i,A) 


ui(A)-C'i(/i,A) 


-8C2(/2, A)^2(/2, A)^^4WA_wpAl + 4 c2(/2, A) 

OiPl, A) 
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Arguments of ^ 



Figure 3: Case /i = 1, ^2 = 1 


After some clean-up we get 

?>Sl{l2, A)/if(A) - 8 A)C2(/2, A)^2(/2, A) + Slih, A)C'i(/i, A)] 

2Siil2, A) - ASfih, A) + 4^2(/2, X)C^Ah, A) + ACHk, X)Sl{h, A) 

+8S^{h, X)C2{l2, X)S2{l2, X)C,{h, A)] fxliX) - Slik, A) = 0. 

The equation satisfied by fi 2 {X) is obtained by interchanging the 
subscripts 1 and 2. 


5.2 Numerical work 

Figures [^1^ and [^display multiplier data for three cases. In all cases 
g = 0 and li = 1. The values of I 2 are: (i) I 2 = 1, (ii) h = 0.89, and 
(hi) I 2 = 2. 

For a range of positive values of a, solutions of the degree four poly¬ 
nomial equations for /ri(cr) and /i 2 (o') are computed. Actual multiplier 
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pairs (/ii,p 2 ) must satisfy the system (4.8), as well as the bounds im¬ 
plied by the square integrability condition (4.11). To eliminate spuri¬ 
ous solutions, the expressions in (4.8) were evaluated, and candidate 
pairs (/ii,/r 2 ) were rejected if either equation had an expression with 
magnitude greater than 10“®. Pairs were also rejected if either candi¬ 
date multiplier had \^j\ > 1, or if the minimum multiplier magnitude 
exceeded 1 / a/3 . 

Each hgure contains two parts, the multiplier arguments and the 
logarithm of the magnitudes. Figure is the case with h = h = 1- In 
this case the two multipliers are equal. By Corollary 4.15 real points 
in the resolvent set can be recognized by real values for both multi¬ 
pliers, except when a lies in a discrete exceptional set. Eigenvalues in 
these sets are possible, as discussed in [5]. 

Figure illustrates the case I 2 = 0.89. When the multipliers are 
not real they will appear in conjugate pairs. Unlike the classical Hill’s 
equation, multipliers may vary in magnitude when they are not real 
valued. The multiplier arguments may exhibit occasional discontinu¬ 
ities. 

Figure illustrates the case I 2 = 2. The multiplier argument dis¬ 
continuities are clearly visible. Notice that the horizontal axis displays 


the predicted periodicity from the proof of Corollary 4.17 

dent. 


IS evi- 
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